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ABSTRACT 

Descriptions of several digital computer programs for electronic 
circuit analyses are given. A listing of a special-purpose (Class '"'C!! 
amplifier design) program is given. Comparisons of several general 
purpose peaneena™ are made. An extensive description of the CADLIB 
is included. The discussion of the CADLIB points out numerous 
modifications to the general purpose programs. These modifications 
include such things as the addition of Minter-plots , more accurate 
algorithms, etc. Also many of the problems encountered in adapting 
the various programs to the NPGS machine (IBM 360 model 67) are 
enumerated. Appendicies are included to provide a reference for coding 
the members of the CADLIB for which no instruction manual was avail- . 
able. An extensive bibliography on applied computer aided design and 
the job control language details of how the CADLIB was created are 
included. In addition two extensive examples using the CADLIB are 


included. 
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CHAPTER I 


INTRODUCTION 


In circuit design a most time-consuming, error-prone, and 
frustrating chore is the repetitious solving of network equations or 
the repetitious measurements taken from a breadboarded circuit as 
some parameter(s) is varied. However, the digital computer thrives 
on simple-minded repetitious work while attaining perfection in the 
writing of equations and near perfection in their solutions. The 
accuracy of the solutions may be even further enhanced by the use of 
double precision calculations which, under the present trend of com- 
puter design, is very nearly as fast as single-precision calculation, 
e.g., the ADD operation in the IBM 360 requires 2.43 microseconds 
in single-precision and 2.45 microseconds in double-precision. 

Another point in regards to accuracy is that no solution is any 
more accurate than its mathematical model irrespective of the number 
of significant digits carried in the computation process. Modelling 
of resistors, coils, and capacitors with the necessary fringing or 
stray effects presents few problems. However, an accurate model 
for an MOS-FET or integrated circuit is prohibitively complicated 
for the slide-rule or calculator approach. But by making use of the 
speed of the computer virtually any model can be effectively used. 

In the following discussions a basic knowledge of the Fortran 


programming language is assumed. There are also numerous 


references to OS/ 360 job control language but the understanding of 
these references is not necessary to make use of this report. In 
particular the discussions on how the various programs were placed 
on a disk is intended primarily for guidance to one who is interested 


in modifying or adding to these programs. 
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CHAPTER@W 


COMPUTER AIDED DESIGN PROGRAMS 


INTRODUCTION 


There are several hundred digital computer programs presently 
in use for circuit design. These programs run the gamut from simple 
specialized programs for a particular circuit configuration to extremely 
complicated packages that have widely varying capabilities such as: 
sensitivity coefficients, worst case, DC, AC, and transient analysis, 
plotted outputs, and recently on-line graphics (light pen and CRT). 
Appendix A contains an example of a specialized program to ''optimize'! 
a Class ''C'"' amplifier feeding a resistive load. As is characteristic 
of these specialized programs the Class '"'C'' program is of limited 
value and the input 1s very restrictive both in content and form. Several 
worthwhile specialized programs in various Fortran languages have 
been listed and documented in JPL's Technical Report 32-979, 
references 18 and 19. Examples of the general purpose type of 


program comprize the remainder of this report. 


GENERAL DISCUSSION 


As pointed out above the capabilities of the available programs 
vary so widely in applicability and scope that any comparison must be 
limited to the few large general-purpose programs. Table 1] lists the 


major characteristics of some of the major programs. The input to 


1] 


these programs is usually ae avea’ and eC eee 1.@., 
not Fortran. Typically the user is required to: 

(1). Number the nodes sequentially from zero or one, 

(2). Either number the branches sequentially or give them a 


unique name, and 


(3). Code the branches on cards using commas as delimiters. 


The limitations of these programs must also be considered. 
All state-variable based programs suffer from two common limitations. 
These are: 
(1). Numerical errors including round-off errors and 
truncation errors, and 
(2). The time constant problem generated by the network 


being analyzed. 


The first limitation is caused by the limited number of digits in a com- 
puter word and can be minimized by the use of better algorithms and 
double-precision arithmetic. The second limitation is the result of 
solving continuous integral equations by quantized numerical methods. 


To obtain accuracy of integration a small step size is needed for even 


By engineer-oriented we mean that a branch is described by its 
node numbers and a positive direction of current flow. 


By free format we mean that specific columns on the card are 
meaningless and delimiters such as commas are used to separate the 
data. 


HEZ. 


om <a ECAP NET- 1 CIRCUS;|SCEPTRE | PREDICT [INASAP 


Init. Cona. 































Sens. Coeff. 
DC Analysis 
AC Analysis 
mrans. Anal. 
Max Element 
Max Nodes 

Stored Model 
Volt. o0Uu. 


@urr. Sou. 


Language FTN/FAP 
Deck Size 10,000 
Formulation St Var 
Integration Trap 


Can be constant 
Can be function of time, table, or equation 


Can be a function of a subroutine 


Table 1. A comparison of the Major Programs. 


le 


the smallest time constant of the circuit. But to obtain the complete 
transient or steady-state solution the integration must be carried out 
to about five of the largest time constants of the circuit. Thus we once 
again hit that perennial ''trade-off'' problem and about the only solution 
at the present is to eliminate as many of the small time constants as 
practical and/or use a variable step-size integration technique to 
minimize the error. 

Another method of circumventing the first limitation is to use 
a topological approach. However, with growing circuit size computer 
run time rapidly becomes excessive due primarily to the requirement 
that all possible trees must be found. This tree-finding problem will 
be treated more thoroughly in a later section dealing with the CALAHAN 
program. At the moment topological methods are limited to networks 


of 25-30 branches and 15-20 nodes. 


ECAP (Electronic’ Circuit Analysis Program) 

ECAP [16] was developed through a joint effort of IBM and 
the Norden division of United Aircraft Corporation and released in 
June 1965 for the IBM 1620 computer. It has since been translated for 
use on the IBM 7040/44, 7090/94, and the System 360 series. ECAP 
with modifications is probably the most widely used general purpose 
CAD (Computer Abb Design) program today. It is one of very few 
(if not the only) general purpose programs to offer both DC and AC 


analysis with an automatic parameter modification feature in both its 
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DC and AC analyses. ECAP will provide, as output, the sensitivity 
coefficients, worst case, and standard deviations all of which are 
noticeably lacking in many other general purpose programs. Another 
positive feature is the switch (not available on most other programs) 
which provides a very flexible piece-wise linear modelling capability 
to approximate non-linear elements. 

A major limitation of ECAP is its lack of a non-linear capability 
except for its extremely good piece-wise approximation capability. 
The major problems in using ECAP are its lack of a model library for 
frequently used components and the problem of data processing, e.g., 
a single node voltage is unobtainable but must be separated from a 


tabulated list of all node voltages. 


NET-1 

NET- 1 [10] was developed at Los Alamos Scientific Labora- 
tory of the University of California at Los Alamos, New Mexico, under 
the auspices of the United States Atomic Energy Commission. NET-1, 
like ECAP, is a general purpose program but suffers most noticeably 
from the lack of an AC capability and the unfortunate (?) choice of 
FAP as a primary programming language; thereby, limiting its users 
to the IBM 7090/94 and Maniac II computers as conversions to other 
machines would be a major undertaking. However, NET-1 does 
have a built-in non-linear model for transistors and diodes but these 
models require some 36 parameters for the transistors and 13 for 


the diode and the model must be entered before it can be used. These 
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parameters are not readily available from manufacturers and must be 
measured from several devices of each type to obtain typical values. 
NET-1 became available in October 1962 and has since been rather 
widely distributed with its users swapping semiconductor character- 


istics rather freely. 


CIRCUS 

The CIRCUS program By & 13] was developed at the Boeing 
Company, Seattle, Washington, primarily for use in studying the 
effects of nuclear radiation on electronic circuits, but its applicability 
to transient analysis of all types led to its use as a general purpose 
program. CIRCUS has the modeling capability (and limitations) of 
NET-1 but was written in Fortran IV (98%) and Assembly (2%) languages, 
and the conversion to the System 360 has been made by Lincoln Labora- 
tories at MIT. CIRCUS is relatively fast and accurate but suffers 
from the same fixed models as NET-1. This can lead to unnecessarily 
long run times due to the small internal time constants of some charge- 


control models. 


SC Dy sadly le, 

Like CIRCUS a desire for a transient analysis of the effects 
of nuclear radiation led to the development of SCEPTRE [11] . It 
was developed by IBM for Kirkland AFB. One of the contract stipula- 
tions was ''100% Fortran'! and this led to a very machine dependent 


program completely defeating the purpose of the Fortran restriction. 
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SCEPTRE has two unique features not found in other programs. 

The first is the ability to scan the loop matrix of the circuit, 
select an efficient partition, generate a Fortran program describing 
the circuit, and then execute the generated program. If desired, this 
Fortran program can be saved and executed at a later time or in steps 
as might be desired when the run time becomes long. Parameter 
modifications make changes only in this generated program so that the 
load time (time spent generating the Fortran program) is encountered 
only once for each circuit. This load time originally required 4-5 
minutes of run time; however, this has been lowered to about 2 minutes 
[21] through the use of more efficient overlays. 

The second unique feature of SCEPTRE is its tremendous 
modelling flexability. SCEPTRE will accept and store any circuit in 
its model library. The connections and their order is specified at 
the time of entering the model into the library and a simple call will 
connect the model anywhere it is desired. The value of the model 
library for ''standard'" or often used circuits is readily apparent. Due 
to its size and its modelling capability SCEPTRE is probably potentially 


the most powerful general-purpose CAD program available today. 


PREDICT 
PREDICT is the forerunner of SCEPTRE but as originally for- 
mulated it has no stored model capability and had a rather simple- 


minded trapozoidal integration algorithm. However, an IBM 360 


lhey 


version has been written and some of the features of SCEPTRE were 


added to PREDICT to form PREDICT- 360. 


NASAP (Network Analysis for System Application Program) 

NASAP [4 & 7] is mentioned primarily because it has had a 
maximization of advance publicity in spite of the fact that it is not 
operational on any machine at the moment. But NASAP is planned to 
have virtually all of the features of most of the available general-purpose 
programs available today. NASAP requires no lengthy matrix inversions, 
thereby eliminating the primary source of round-off and truncation 
errors. NASAP utilizes a dichotomous approach with a topological 
analysis of many small circuits and preliminary testing indicates that 
the computer run time is reduced by factors as large as 10° from the 


state variable based programs. 
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CHAPTER Ill 


CADLIB 


The CADLIB is a computer library of computer programs 
designed to aid the designer in the analysis of various circuits. The 
CADLIB is presently maintained in executable form ona direct access 
disk named J@BO003 in the computer center of the Naval Postgraduate 
School. Support for the CADLIB is furnished in the form of a sequential 
data set named MODEL (also on J@BO03) for use as a model library for 
one of the members of CADLIB. The various members of CADLIB are 


executable members of a partitioned data set and they consist of the 


following: 
ECAP slightly modified by J. G. Sketoe to furnish 
printer plots fm . 
CIRCUS OS/ 360 version distributed by MIT except 


that the default model unit 11 capability 
was inserted. 


CALAHAN slightly modified by J. G. Sketoe to furnish 
automatic printer plots both of freq and 
transient responses. 


POTTLE slightly modified by J. G. Sketoe to furnish 
automatic printer-plots of the freq response. 
Also the matrix inversion has been changed 
to double-precision arithmetic. 


rPOorror the POTTLE program with a free-field format 
capability as designed by J. G. Sketoe 
(See Appendix D). 


MODEL the CIRCUS modeling program as described in 
reference 13. Operates on data set MODEL. 


bas 


The CADLIB data set was constructed using a SYSLM@D of 
the following form: 


//SYSLM@D DD DSNAME=CADLIB, UNIT=2311, DISP=(NEW,KEEP), xX” 


Ihe: V@LUME=SER=JQB003, xX 
ei SPACE=(CYL,(25,2,2)), x 
bai LABEL=EXPDT=99350 

And the individual members were added using 
//SYSLM@D DD DSNAME=CADLIB, UNIT=2311, DISP=QOLD, xX 


V@LUME=SER=J@BOO3 
Execution of the various members requires different sequences 
of control cards. These sequences will be given below as each member 


is discussed. 


meoAP 





The first program obtained for the CADLIB was IBM's Elec- 
tronic Circuit Analysis Program. The 360 version requires approx- 
imately 290K bytes of core when used without overlays as in CADLIB. 
ECAP was obtained from Share, a facility of IBM's, in the E-level 
compiler edition. As the OS/ 360 version was ordered only the 360 
manual f20] was included. However, this manual only pointed out the 
differences in ECAP-1620 and ECAP-360. A manual was then ordered 
for the 1620 version (6) and two months later it arrived. In the mean- 


time by using ECAP problems published in various journals it was 


ale 
~~ 


“This X as in all of the following control cards must be in column 72. 
Also the first character of the following card (except for the double 
slashes) must be in column 16. 
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discovered that the received version simply would not work properly. 
The major problems and their solutions were: 

(1). The release 9 G-level 'in-house'' compiler failed to con- 
struct the proper instructions for overlay. At that time 
the computer center policy was to run two-256K byte 
machines rather than one 512K byte machine so that 
various hardware testing could be done concurrently with 
computer service. Use of the older but slower H-level 
compiler cured this trouble. (This fault of the G-compiler 
has been cured in later releases). 

(2). The next problem discovered was that the language 
processor portion of ECAP did not function correctly 
when compiled on either the G- or H-level compiler. The 
cause of this trouble lay at the fundamental differences 
in the E-level Fortran compiler and the G/H level com- 
pilers. As the language processor converts an engineer- 
oriented circuit description into the required computer 
description, literal data (alphanumeric characters) is 
manipulated vice purely binary numbers (floating point 
numbers are a special coding of binary). This manipula- 
tion is done by using a unique combination of bits to 
represent a particular character. When an arithmetic IF 
statement is used to decode this data, the data is treated 


as numbers and not as literal characters. For example, 


2) 


(3). 


the Fortran statement 

IF (M-N) MO 15; 10 
compiles into an Behm etic subtraction of N from M, 
examine the resultant, then branch to statement 10, 15, 
or 10 according to whether the resultant is negative, zero, 
or positive. The contents of a computer word or cellis 
considered positive or negative according to whether the 
first bit is a zero or one, respectfully. In the IBM 7090/94 
machines the literal data codes were mixed as positive 
and negative numbers with the special characters (,. etc.) 
intermingled within. This mixture was observed inv thie 
E-level machines while in later machines this coding was 


changed such that alphanumerics are negative (and 


sequential) and special characters are positive. The 


solution to this problem was quite obvious once the source 


re Fj 


of the aS was known. All literal data arithmetic IF's 
were ee deed with logical IF's which do their comparisons 
- a pH Ser - bit basis rather than an arithmetic answer. 
senine of — two problems finally got our copy of ECAP 
at least partially Wonidine. The next problem encountered 
me that the automatic parameter FaTIvaSE feature did 

not function properly. A cry for help to IBM was answered 
with a list of aaa and omissions in the originally dis- 
tributed erin of ECAP. The appropriate entries cured 


this problem. 
a2 


(4). Next an addressing error was encountered when trying 
to obtain a frequency response solution with an overlaid 
program. As the various control sections in overlays 
complicate debugging, an attempt was made to find the 
error while using all 512K bytes of core without overlays. 
The addressing error disappeared. Reinsertion of the 
overlay brought the bug back. The only conclusion was 
that the bug was in the machine constructed overlay. 
Next run time comparisons were made on overlay vs. 
non-overlay versions. The results were that overlaid 
programs required 2-4 times as much run time as non- 
overlay. By this time all 512K was being used for almost 
all runs and the decision to use ECAP without overlay in 
the CADLIB was made. The bug was brought to the attention 


of the computer center personnel and left there. 


Presently all known bugs have been worked out of the NPGS 
version of ECAP in addition to the addition of a printer plotting addition 
as explained in reference 17. As User's Manuals are in short supply 
an abbreviated manual has been written and is published as reference 17. 

The required control cards to execute the CADLIB version of 


ECAP are shown in figure 1. 


ZS 


//ECAP J@B (EAC6,01FT),' NAME ',MSGLEVEL=1 
//J®@BLIB DD DSNAME=CADLIB, UNIT=2311, DISP=@LD, 
le V@LUME=SER=J®B003 

fea EXEC PGM ECAP 

//FTO6F001 DD SYSQUT=A 

//FTOS5F001 DD * 


Batched ECAP coded input decks 





Figure 1. Control Cards to Execute ECAP. 


Basically ECAP is formulated to fill the nodal conductance 
matrix as the input branches (B-cards) are readin. The various 
command cards then set am appropriate flags to indicate which 
algorithms are executed to obtain the desired outputs. This conductance 
matrix is available as an output request under the MISCELLANEOUS 
request. 

One additional note worthy of comment about ECAP is that it 
satity interprets the first two characters of any word, e.g., PRis the 
same as PRINT, ‘e is the same as V@LTAGES, BE is the same as 


BETA, etc. 


CIRCUS 

A. copy of the CIRCUS-360 program was obtained from Lincoln 
Laboratory, MIT, Lexington, Massachusetts. As with ECAP the 
received version would not Work although MIT had exactly the same 
type machine as NPGS, i.e., a model 65/67 IBM 360. The MIT des- 


cription of the distributed CIRCUS tape is given in Appendix B. The 


various problems encountered with our copy were: 


ee 


(2). 


(Se 


The CIRCUS tape, as received, had an H-level compiler 
output version on file # 1 (Beene B) which should have 
loaded directly from the tape. However, a tape error 
had accu eel and the tape had to be recompiled. 

At MIT there exists three resident routines that access 
the day, year, hour, minute, and second. These 
routines were used by CIRCUS-360. As these routines 
did not exist at NPGS the initial solution was to provide 
dummy routines to get CIRCUS operating. While this 
solution worked it eliminated the variable time limitjon 
CIRCUS as explained in the User's Manual 112] . The 
author then proceeded to learn enough assembly language 
to write a set of routines to accomplish this. These 
routines are now added to our version of CIRCUS as well 
as entered in the resident library of our computer center. 
While translating CIRCUS MIT evidently inserted a de- 
bugging print-out at each integration step. As this led 
to excessive and useless print-out of the integration step 
time and size, both in time and paper this routine was 


slightly modified to eliminate this waste. 


“A bEND card had become BEND where "'b'' implies a 12-2-9 punch. 
“Ao Character implies a 12-2 puncn. 


ZO 


(4). 


(5) 


No overlay table was included in CIRCUS- 360 as it will 
fit into a 256K byte machine if the load point is below 
BO000 (See Appendix B). However, the load point at NPGS 
is B820¢ with just the "in-house'! G-compiler thereby 
necessitating an overlay. This overlay was constructed 
and is shown in figure 2. This overlay was used but was 
soon dropped for the following reasons: 

(a). Substantial decrease in run time. 

(b). Prior experience with ECAP indicated 
GErOrs Mimeht Oecur. 

(c). Using the 256K byte machine only 158K 
was available with HASP. CIRCUS even 
with this overlay required 153K leaving 
almost no room for increased compiler 
size (which came later). 

(d). By going to 512K we were able to increase 
our maximum circuit size to circuits of 


over 70 equations (Appendix B). 


The MIT description lists two known bugs (Appendix B) 
in their program. The second one of these failed to show 
up in our recompiled version of CIRCUS. This led toa 
successful attempt to correct the tape error mentioned 


earlier so that this bug could be checked out. As indicated 
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Root 
Segment 


CHANGE 
CREATE 
DICT 
ELFIND 
FETCH 
GENRAT 
INCIDN 
JFIND 
nF IND 
LINKI 
LINK2 
MAPPER 
PLSTOR 
REGIT 
SETUP 
SHIFT 
SMODEL 
STORV 


MAIN 
CLOCK 
ENDJOB 
EQUATI 
KLOCK 
KLOCKI 
MATOUT 
MONKEY 
OUTMAT 
TIGER 
COMPR 


BEAR 
ELPHNT 
HORSE 
INDUCT 
LINK4 
TMPLY 
TRIPLE 


Figure 2. 


MAINZ 
DOUBLE 
MPLY2 
rae OE: 
SOURCE 
SVALUE 
TABLIZ 


CHARET 
EXPINT 
GAMT 
INTRPT 
IS ee db, 
LINK5 
LISTVL 
NODEV 
Oe ar 
PHOTO 
PLOT 
PRINT 
Pies 
SPCAL 
SPeArD 
VMPLY 
XXMIN 


oa 


EAS 


UNPAC 

EXTRCT 

GET 

PUT 

BCDEP'T 

SETLMT 

BCDINT 

TIMHMS 

TIMHR 

TIMDAT 
FINDER 
LINK6 
Jo) of Oud i Dyas 
CLEAR 
Peres! 
PUPP L 
FINDY 
SCALEY 
FINDX 
SCALEX 
SPTPLT 
PLTOUT 


CIRCUS- 360 Overlay Requiring 
approximately 153K bytes. 


Overlay 
Alpha 


Overlay 


Beta 


in Appendix B the MIT version failed to recognize a 
card of the form 
USE M@DEL LIBRARY, n 

where nis integer denoting the logical unit of the model 
library. The NPGS version did recognize the card 
correctly. Naturally the NPGS version was placed in 
CAD Le 

(6). An error of omission by MIT but not mentioned by them 
in their write-up was that the default to unit 11 (page 25 of 
User's Guide [12] ) was omitted. This has been comme 


in the NPGS version. 


Presently CIRCUS- 360 is in CADLIB and is working properly. 
As CIRCUS uses several different logical units a cataloged procedure 


has been established for CIRCUS in order to simplify the job control 


card requirements. The standard execution of CIRCUS using logical . 


unit 11 for the model library, when needed, is shown in figure 3. 


//CIRCUS J@B (EAC6,01FT),' NAME ',MSGLEVEL=1 
//J@BLIB DD DSNAME=CADLIB, UNIT=2311,DISP=@LD, 
~ VQ@LUME=SER=JQ@B003 
EXEC CIRCUS 
//GO.SYSIN DD * 


Batched CIRCUS coded input decks 





Figure 3. Standard CIRCUS Job Control Language. 
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If logical assignments other than or in addition to the default 
assignments are desired (such as use of 'SAVE TAPE"), these units 
must be defined in a manner similar to that shown in figure 4 where 
unit 15 is defined for use asa "SAVE TAPE! or as a model library tape. 

A listing of the actual cataloged procedure to execute CIRCUS 
is given in figure 5 for reference althouth it should not be needed 


except to perhaps extend or replace it. 


//GIRCUS J@B(EAC6,01FT),' NAME ', MSGLEVEL=1 

//IJGBLIB DD DSNAME=CADLIB, UNIT=2311, DISP=Q@LD, 

// V@LUME=SER=J@B003 

// EXEC CIRCUS 

//G@.FT15F001 DD DSNAME=SAVE, UNIT=2400, DISP=(NEW,KEEP), X 
V@LUME=SER=XXX, LABEL=YYY, xX 


DCB=(BLKSIZE=1204,RECFM=V) 
//G@.SYSIN DD +* 


Batched CIRCUS coded input decks 





Figure 4. Defining Logical Unit 15 for use by CIRCUS 


MEMBER NAME CIRCUS 


//GO EXEC  PGM=CIRCUS 

//FTOIF001 DD DSNAME=SYS1.SYSUTI1, DISP=@LD, x 
a DCB=(BLKSIZE=1204, RECFM=V 

//F TO3F001 DD DSNAME=SYS1.SYSUT2, DISP=Q@LD, x 
// DCB=(BLKSIZE 1204, RECFM=V) 

//FT12F001 DD DSNAME=SYS1.SYSUT3, DISP=@LD, xX 
‘A DCB=(BLKSIZE=1204, RECFM=V) 

//FTUFO00l] DD DSNAME=MQDEL, UNIT=2311, DISP=Q@LD, xX 
had V@LUME=SER=J@B003 


//FTO6F001 DD SYS@UT=A 
//FTO5F001 DD DDNAME=SYSIN 


Figure 5. Cataloged Procedure "CIRCUS" 
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The formulation of CIRCUS is based in the state varible approach 
as advanced by Bryant [22 | and extended by Wilson and Massena (i235) 
A detailed, although not rigorous, explanation of the mathematical 


theory of CIRCUS is given in the Program Manual fis ' 


MODEL 

The MODEL member (alias UPDATE) of CADLIB is a program 
to list, modify, generate, etc. the model library as described in 
Chapter 4 of the CIRCUS Program Manual [13] . The CIRCUS tape 
as received from MIT had a compiled version of the MODEL program 
as file #2 (see Appendix B). This file was used as input to the linkage 
editor to form member MODEL which operates on the actual CIRCUS 
model library. The control cards required for execution of the member 
MODEL is shown in figure 6. The model library is described by the 
FT11F001] card in figure 6 which, as shown, describes a date-protected 
sequential data set named MODEL which resides on the same JOB003 
disk as the CADLIB. If one wishes to use the program or member 
MODEL to operate on some other data set (i.e., magnetic tape) only 


this F T11 F001 must be changed. 


au) 


//M@DEL J@B (EAC6,01FT),' NAME ', MSGLEVEL=1 
//IJ@BLIB DD DSNAME=CADLIB, UNIT=2311, DISP=Q@LD, 
V@LUME=SER=JQ@B003 
PGM=M@QDEL or UPDATE 
//FT11F001 DD DSNAME=M@Q@DEL, UNIT=2311,DISP=@LD, xX 
a V@LUME=SER=J@B003 
//FT12F001 DD DSNAME=SYS1.SYSUT1, DISP=@LD 


//FETO6F001 DD SYSQ@UT=A 
//FTOS5FO0Ol1 DD * 


Model input as per Chapter 4, CIRCUS 
Program Manual [13] 





Figure 6. Use of the Model Program. 


BOTTLE 

The copy of C. Pottle's Linear Network Analysis Program as 
received required only changes in the logical unit assignments to make 
it compatible with the NPGS machine. To make the use of this program 
feasible from the CADLIB all outputs were assigned to unit # 6 (printer) 
and the input was assigned to unit #5 (card reader). However, flexibility 
was added to the program by making the unit designations variables 
in the individual read and was it, aisha ieee and then assigning these 
variables in a central location in the main program. These variables 
are then transferred between various subroutines via labeled common. 
This provides the user with the facility of writing out some data 
separately from other data. The variable names and their assignments 


are Shown in figure 7. 
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Use 

INTAPE User's input data 

IOUTP Copy of data set & errors 

Ss TP | State equations 

TRANTP Transfer functions, poles, & ae 


FREQTP Frequency response 





Figure 7. Default POTTLE I/O Assignments. 


The control cards required to execute the POTTLE program 


from the CADLIB is shown in figure 8. 


//P@TTLE J@B (EAC6,01FT),'' NAME '!,MSGLEVEL=1: 
//J@BLIB DD DSNAME=CADLIB, UNIT=2311,DISP=Q@LD, 
V@LUME=SER=JQ@B003 
PGM=PQOTTLE 
//FTO6F001 DD SYS@UT=A 
//ETO5F001 DD * 


POTTLE coded input (See Appendix C) 





Figure 8. Control Cards to Execute POTTLE. 


The procedure |14] implemented by Mr. Pottle is based ona 
recursive algorithm developed by Derisoglu C6 | which extends Bryant's 
algorithm to nonlinear active elements. The state equations are formula- 
ted;.then the frequency response is determined via the Souriau-F rame 


= 
algorithm to obtain the matrix (sI-A) . The time responses are 
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obtained by factoring the characteristic polynomial and performing 
partial fraction expansions of each element of the transfer function 
matrix or by direct integration of the set of first-order equations. 

The instructions for coding the POTTLE program are included 
as Appendix C. 

Changes to the POTTLE program were limited to (1) providing 
extended accuracy through the use of a double-precision vice the 
original single-precision matrix inversion routine and (2), the addition 
of an automatic plot of the time response and the linear frequency 


response. 


POTTOL 

The frustrations encountered in using the POTTLE program 
pointed out the desirability of a translator so that a free-format input 
could be used. The author has designed such a translator and combined 
it with the POTTLE program described above to form the POTTOL 
program. While developing the Translator several other convenience 
features were included. 

While the use of the POTTOL input format is relatively simple, 
a complete description of the Translator is not so simple. The complete 
description of the Translator input coding is given in Appendix D. The 
control cards to implement POTTOL are shown in figure 9 and are 


very similar to those used for ECAP and POTTLE. 


3) 5, 


//P@TT@L JIOB (EAC6,01FT),' NAME ',MSGLEVEL=1 
//J@GBLIB DD DSNAME=CADLIB, UNIT=2311,DISP=Q@LD 
i V@LUME=SER=J@B003 

; rx eC PGM=PQOTTOL 

//FTO6F001 DD SYSQ@UT=A 

//FTO5F001 DD oe 
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Batched POTTOL coded input (See Appendix D) 





Figure 9. Control Cards to Execute POTTOL. 


CALAHAN 

Probably the most accurate program in the CADLIB is the 
CALAHAN onl which was developed by D. A. Calahan while at the 
University of Illinois. As received, the program was coded in IBM 
7090 Fortran II, The problems encountered with the program were 
generated primarily by a complete lack of knowledge of the intricacies 
ereF ortramell. But = differences were not as great as feared. | The 
names of the Fortran library routines were slightly different (e.g., 
COSF vs. COS) and the passing of arguments through subroutines eet 
be handled differently causing the computer to get lost when passed 
in the same manner in FortranIV. The major problem encountered 
in implementing and wettnie the program was the intepretation of the 
instruction manual. siepieoneltm E is a set of instructions which (hopefully) 
will clearify the coding of the input to the CALAHAN program. 

The CALAHAN program is a topologically based program which 


is the primary source of its accuracy. This approach evaluates the 
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determinants and cofactors of the node admittance matrix by finding 

the related n-tree products of the corresponding linear graph. Basically 
the matrix solution for the node voltages is transformed by the Binet- 
Cauchy theorem 115 ] and the Maxwell's formula to a problem of 
determining all tree and associated n-tree products. This tree-finding 
problem quickly replaces the integration time constant problem of the 
state-variable program as a computer run time limitation. Regardless 
of the tree-finding scheme used the possible combinations of branches 
increase factorially (about the most rapidly increasing function known) 
and the number of trees increase exponentially (about the second most 
rapidly increasing function known). For example, the circuit of figure 
10 using the hybrid-pi transistor model results ina circuit containing 
eleven nodes and twenty-eight branches which forms 456,172 trees! 
While all topologically based programs have this tree-finding limitation, 
they do not suffer from the round-off error so prevalent in the state- 


variable formulations. 


OV ae 


oe =a 


Figure 10. Wide-band R-C Coupled Amplifier. 


SS 


The execution of the CALAHAN program from the CADLIB is 


shown in figure 1] and again the similarity to the execution of the 


ECAP and POTTLE should be noted. 


//CALAHAN JOB (EAC6,01FT,' NAME ',MSGLEVEL=1 
//J@BLIB DD DSNAME=CADLIB, UNIT=2311, DISP=@LD, 


fe! 
i 


V@LUME-SER=J@B003 
EXEC PGM=CALAHAN 


//FTO6F001 DD SYSQ@UT=A 
//FTO5F001 DD * 


ere, 


Batched CALAHAN coded input (See Appendix E) 





Figure 11. Control Cards to Execute CALAHAN. 


The extensions and additions made to the CALAHAN program 


(1) An automatic printer-plot of the requested output for 
both frequency response calculations and time response 
Cai cumenul Ons, 

(2) A double-precision root-finding algorithm which 
usually works better than the original ''Method of 
Mueller" algorithm has been added. The roots are 
found by both methods, two polynomials are then con- 
structed from the two sets of roots, the magnitude of 
the coefficients are compared to the original. The set 
of coefficients with the least error determines the 


selected set of roots. 
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Where 








V. 
ZA Sia SES) = _ 
hal 
And 
Noa 
T(s)s= @E Cs) = 7 
in 


Considering now the case of non-ideal inductors and capacitors 
we can define a new function or ''s'' in which the impedances/admittances 


become (r, + Ls) / (Se + Cs) where r_ and 8, are the dissipative 


Mey 


elements of the inductors and capacitors, respectively. Now defining 


_, ee _ Jg 
o> TL = 


We have 


25) = 50S) aud ¥gés).= C.Qscuiglh 


At this point we could depart from the more rigoreous approach 
and make the assumption that 3, > oe ss ; thereby, making life 
very simple by drawing the conclusion that our function with non-ideal 
elements is formed by the substitution of (s + § ) for s in the original 
function. In order to test the validity of such an assumption the 
dissipation of capacitors and coils were measured using a General 
Radio Type 1650 impedance bridge. These components were chosen 


at random with the only criteria being that a "usage range" be obtained. 


oye 


The results of these measurements are tabulated along with calcula- 


tions of 3 and % intable 2. As is apparent from the table the 


Type Nomial Measured D( ) Cal. Comments 

Cap Electro 4 ufd 4 ufd | 270 lice Ignored 
lyre see text 

Oil 2 wid Ns aay 0.0245 

Paper . 001d SO0093a. 1.2 0.0113 Minimum 

Mica 0047 .00446 4.7 0.0295 

Mica 200 pi W2 pi 4.0 0. 026 

Cermic 250 pf 258 pf 120 0.0628 Maximum 


oorun 7 ieeoh 
Minimum 
RF - Choke 
Air- Dux* 
Air - Dux* ; : Maximum 


Locally 
wound 


* Air -Dux is a trade name of a type of coil used for 
transmitter-to-antenna lead coupling 





Table 2. Coil and Capacitor Measurements 


numerical value of § or S. did not vary excessively over the ‘usage 
range" except for the single 4 ufd electrolytic capacitor. As electrolytic 


capacitors have relatively little use where the pole and zero positions 
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are critical, this particular capacitor was subsequently ignored. 
However, the approximation that 2,=, must certainly be held invalid 
or at least suspect as these parameters differ by factors of 100. 

By utilizing the digital computer, the analysis of a particular 
circuit can be reduced to the coding of a few cards. Utilization of 
the CALAHAN digital computer program in accordance with the 
instructions in Appendix E provides the desired quantities with a 
minimum of druggery. 

The two circuits used for this problem were: 

(1) The Pi-Section of figure 12 which is often used to couple - 

a transmitter to an antenna, and 


(2) The Chebyshev filter of figure 13. 


Figure 12, Pi-Section Coupler 


In. /h 2h) 


|o.334 let Taos Vout 


Figure 13. Chebyshev Filter 
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.505 + jl.321 2194+ jl. 


525 + jl. 313 .2368 + jl. 
550 + jl. 303 2585 + jl. 


. 750 + j1.200 4342 + jl. 





Table 3. Pi-Section Pole Movement (r, = 0, vary Je ) 


Since it was suspected that the capacitors would be substantially 
“more ideal'' than the inductors in their effects on the circuit the first 
computations were to vary Be over the ''typical range'' andydetermine 
the functions' pole positions of the Pi-Section. These are shown in- 
Table 3. Next similar computations for the ''typical!' Hi range 


were made. These results are shown in Table 4. 


-1.0 + j 1.4142 


=i S68 , -4.414 


-1.228, = Oem c 





Table 4. Pi-Section Poles vs. 


As can be seen from these calculations the effects of the typical 
inductor on the movement of the poles is some 100 times greater than 


the typical capacitor. This correlates nicely with the measured § ¢§ 


4] 


The computations for varying De in the Chebyshev filter 
of figure 13 were carried out but the effects were again on the order 


of 2% (See table 5). 


-0.005 + j1.514 





Table 5. Chebyshev Filter (Poles and Zeros vs Me 
© 
Cl C2 R 
V ; Vout 
nm if if In. Ou 
© 
/ 
sy ast 1 
52+ (RCI+14+ Reo )S + RL +55) 


Figure 14. Pi-Section with Non-Ideal Elements. 


The transfer function of figure 14 was calculated and the poles 
determined for various values of RCl and RC2. The particular values 


were 10, 100, and 1000 ohms which corresponded toa oF ro) ca © 
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0.01, and 0.001, respectfully; thereby, bracketing the measured 
typical values. With each of these values RL was varied over a very 
wide range of 10°! to 100 and pole positions calculated. These results 
are tabulated in table 6 and graphed in figure 15. 

As indicated by the graphs the pole positions certainly do not 
move straight in the negative real direction. In fact, there is quite 
a Severe vertical movement. However, notice that the poles do come 
back down to level with the ideal at RL ~ 2 which is in the midrange 
of the measured é. ; thereby, in some sense justifying the ''direct 
left shift'' as a ball park design criteria. 

The effects of the capacitor admittance is to move the pole 
position slightly to the left and downward while the inductive resistance 
moved the pole up until RL = L (numerically) and then the effect be- 
comes quite the opposite. The major point to be gleaned from this 
problem is that the usually neglected — resistance (in filter 
design) can be of such a magnitude to cloud the analysis (with ideal 
elements) in doubt. The CADLIB provides a tool to perform the 
gucmatical druggery so that these parasitic elements can be included 


im thie Oraginal design. 
ON COMPENSATIONS IN INTEGRATED CIRCUIT DESIGN 


Integrated circuits has become one of the most important areas 


‘in space and military electronics. As a consequence, ingegrated 


circuit design has become one of the major branches of network 
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RL RC1=RC2=10 RC 1=RC 2=100 RC 1=RC 2=1000 


10 + + £5005 + 
10-6 + + 5005 + 
nee + + .5005 + 
0.1 + + £5505 + 
0.2 + + 6005 + 
0.4 as is HOS eecae 
0.5 + + 7505 + 
0.7 + + 8505 + 
0.9 + + -.9505 + 
1.0 + + -1.000 + 
laa + -1.055 + j -1.050 + ; 
an? + -1.105 + j -1.100 + j 
rege + “1.155 + | “1.150 + 3 
1.4 cf -1.205 + j -1.200 + ; 
1.5 + -1.255 + j -1.250 + j 
vee + -1.305 + j -1,300.+ ; 
ee + ~1.355 + j ~1.350 + j 
1.8 nm -1.405 + j -1.400 + j 
1.9 + -1.455 + 3 -1.450 + 
2.0 + -1.505 + j -1.500 + 
2 + -1.605 + j -1.600 + 
2.4 + -1.705 + | -1.700 + 
2.6 + -1.805 + j -1.800 + 
2.8 + -1.905 + j -1.900 + 
3.0 + -2.005 + ; -2.000 + 
Be -2.555; -1.659 5 ore le 
10.0 é ; ~9 8873; -1.123 mm=9..887 20-1. 
50.0 -49.98; -1.030 -49.98: -1. 
100. -99.99; -1.020 599995 ~1Le 





Table 6. Non-Ideal Pi-Section Poles (Figure 14). 
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synthesis. But solid state devices are extremely vulnerable to the 
effects of radiation, since the operation of most of these devices is 
determined by the density and gradients of the minority carriers. 
This minority carrier density can be drastically upset by radiation 
exciting electrons from the valence band to the conduction band of the 
material. Whena transistor ina circuit such as figure 16 is exposed 
to a high energy radiation pulse the resultant collector current wave- 
form is similar to that given in figure 17. 

The integrated circuit transistor operating in the active region 
can be considered basically as one forward-biased and one reversed- 
biased diode connected back-to-back. When irradiated electron-hole 
pairs are generated throughout the semiconductor body including the 
substrate (e.g., silicon). In the transition region near a reverse- 
biased p-n junction the. generated minority carriers are swept across 
the junction by the electric field within a few nanoseconds. Carriers 
generated outside of the transition region will raise the minority 
carrier concentration gradient near the junction generating a diffusion 
current. These two components of current constitute the primary 
photocurrent in a p-n junction and; thereby, a transistor. The general 
shape of the primary photocurrent is very much dependent on the 
minority carrier lifetime as indicated in figure 18. In an npn transistor 
the generated collector region minority holes flow into the base region. 


This hole current is amplified by the current gain of the transistor 
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figure 16, A simple Transistor Circuit 


ovose 








Radiation Pulse 


time 


Secondary Current 


Primary 


Figure 17. Transistor Photocurrents (induced) 
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Relative Primary 
Photocurrent 





©, | 0.2 


Ficure 18. 


Notes Numbers are minority-carrier 
lifetimes in nanoseconds 


OO 
100 
30 


0.3 o4 0.5 0.6 0.7 


time (microseconds) 


Relative Shaves of the Diffusion 
Component of Primary Photocurrent. 
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circuit generating a secondary photocurrent such as that shown in 
faxomere 17. 

Radiation-induced photocurrents can cause errors and failures 
in integrated circuits. If the dosage is large enough the semiconductor 
material can suffer permanent displacements effects in addition to 
the ionization effects which are usually temporary i en] . These 
effects as well as the detailed explanation of the theory and relations 
between the radiation dosage (electron, gamma, etc.) and the resultant 
photocurrent magnitude and shape we leave to those specializing in 
such matters. Of major concern here is the fact that we have 
"parasitic'’ photocurrents whose effects must be at least diminished 
if not eliminated. One method is circumvention, which in its most 
common form amounts to removing power from critical circuits 
during the radiation pulse. However, in many applications (e.g., 


counters, guidance mechanisms) continuous operation is mandatory. 


2. 


*) Yor 





Pupe@@we 19. A Simple Liopic Cate. 


49 


Dielectric (e.g., silicon dioxide) isolation is a very effective 
means of removing many of the "evils'' of radiation. The simple 


logic gate shown in figure 12 is used to explain the effects of dielectric 


Figure 20. Logic Gate with Parasitic Elements 


isolation. If this circuit is built on a silicon chip by a diffused process 
with p-n junction isolation (i.e., isolation between the components 

and the substrate is obtained by reverse-biasing the p-n junction 
boundaries), the parasitic diodes and transistor shown in figure 20 
will generate substantial photocurrents making analysis extremely 
difficult. These photocurrents will vary widely and some will be 

quite large. Also the parasitic transistor forms an SCR and problems 


of latch-up are encountered. By use of thin film resistors and di- 


2) 9) 


electric isolation that is relative insensitive to radiation (wide for- 
bidden region) these parasitic elements and their leakage currents 
can be eliminated and the integrated circuit can be analysed very 
much the same as discrete circuits with, of course, the appropriate 


parameters. The resultant equivalent circuit is again figure 19. 


Vie = 1Ov 








} 
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4 ppd 
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Figure 21. Photocurrent Generator and Compensation 
Diode (dotted). 


The major problem is now the transistor photocurrent, eae 
as shown in figure 21. This photocurrent can be compensated for by 
adding a diode as shown (dotted). If thediode junction is made to match 


the collector-base junction, the photocurrent, i , of the compensating 


ppd 
diode will be the same as _ so that no secondary photocurrent will 


result. The disturbance due to the transient radiation will then be 


kept to a minimum. 
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iF re “sy No Compensation 


Fartial Compensation 


| Primary 
FhotoCurrent 


Figure 22. Compensation Definitions. 


By defining the component compensation as 


) 
% compensation = SS SS ( 100% ) 
N pp 


where the terms are defined in figure 22 it is possible to have some 
measure of the effectiveness of a compensation scheme at eliminating 
the effects of the secondary photocurrents. In the circuit of figure 2] 
itis easily seen that at least some compensation will result if any 
compensating diode is used regardless of whether or not the collector- 
base junction is identical to the diode junction. Also since the photo- 
current of a p-n junction is directly proportional to the junction area 

[ 25] the ratios of Bs and "ond will be the same as the ratios of their 


areas with the amount of compensation being proportional to this 


ratio. 


Sj 


In the circuit of figure 21 the amount of current compensation 
will be determined by how closely the junctions can be made to match. 
However, even with 100% compensation this configuration suffers 
from the voltage drop of — through the collector load which, in itself, 
may very well be intolerable. Also as pointed out by Grannemann, 
et. al., [26] the amount of compensation provided is very sensitive 
to a close match of the compensating diode and the transistor collector- 
base junction. This sensitivity is due primarily to the amplification 
of any mismatch by the gain of the stage. Utilizing the CIRCUS 
computer program {12 & 13} an analysis of figure 21 was made 
using the parameters of the 1N903 and 2N404 as '"'typical'' while the 
compensation diodes were special formulations to match (or mismatch) 
the collector junction of the 2N404. This analysis showed that even 
a primary photocurrent of 0.3 milliamperes peak was intolerable in 
that the output which was 4.23 volts with zero input dropped to about 
zero volts and stayed there for about four microseconds when subjected 
toa i0/° rad per second 200 nanosecond radiation pulse such as shown 
im iaoure IY. 

By attaching an additional identical compensation diode across 
the load resistor as shown in figure 23 we were able to achieve sub- 
stantially better compensation in that are remained more nearly 
constant. The sensitivity of the output to a mismatch of this load 
diode is much less than the base as there is no gain following it. 


But as noted above imperfect matching of the base diode to the transistor 


3. 


Figure 23. Base and Load Diode Compensation 


is magnified. Next over-compensating with the load diode was tried. 
This led to an "increasing" output voltage (with zero input) during a 
pulse. This factor was immediately seized on as a means to decrease 
the sensitivity of the circuit to the amplified mismatch of the base 
diode. The amplified mismatch lowers \ ee while the overcompensated 
load diode raised a This over-compensation can be effectively 
accomplished by increasing junction area of the load diode. By 
utilizing a factor of 2 and our values a compensation of 80% oe 
attainable for a base diode tolerance of 20%. 

Next the circuit of figure 24 was attacked. Initially only base 
diodes on Tl and T2 were utilized and theoretical 100% compensation 
was attainable for each transistor if the diodes and transistors were 


exact matches. However, sensitivity to mismatch worsened as Tl 
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“igure 25. Comoensations for the Gate Circuit. 
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mismatch is amplified by T2 when T2 is turned on by its own mis- 
match. Asa tolerance of less than 1-5% is not feasible in integrated 
circuits with the present technology, this approach was abandoned. 

As T2 forms part of the emitter circuit of Tl attempts to com- 
pensate solely T2 met with very little success. By attaching base 
diodes tolbotH Tl and T2 along with load diodes across the 1.75 K and 
the 6.0 Kohm resistors as shown in figure 25 adequate compensation 
was effected. As mismatch at Tl would cause a secondary photo- 
current into the base of T2 the base diode of T2 should be slightly 
greater in junction area than Tl base diode to provide maximum area 
matching tolerance. With the factor of two between base and load 
diode junction areas a compensation of 80% was attainable with areas 
tolerance of 20%. 

A theoretical development of an effective and economical 
means of radiation compensation in integrated circuits (also discrete 
but not so economical) has been demonstrated. An iterative transient 
analysis is required, thereby, virtually demanding the use ofa 
computerized analysis. The costs associated with constructing, 
testing, and evaluation actual devices without the use of such an 
analysis would be prohibitive while with the analysis at least a reason- 


able first cut can be made. 
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APPENDIX, A 


CLASS "Ce AMPLIFIER 


The following program selects appropriate parameter values 
to optimize a class ''C'' circuit. As rather complete comments are 
contained in the listing a detailed explanation is not included at this 
point, but two points should be brought out. The optimizing function 
1s calculated in subroutine OPT and, therefore, it is easily modified 
to select or weigh some other parameters. Similiarly KCURR 
describes a tube model that can be quite easily changed to some 


other model (e.g., a transistor model). 
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MAP: 


File 1 


File 2 


File 3 


File 4-12 


File 13 


File 14 


FORMAT: 


APPENDIX B 


FILEMAP AND FORMAT 


CIRCUS DISTRIBUTION TAPE 


CIRCUS OBJECT DECK: From FORTRAN H OPT 
O COMPILER and release 9 assembler. 4000- 
5000 cards. May be loaded directly from tape. 


CIRCUS MODEL LIBRARY UPDATE PROGRAM 
OBJECT DECK (complete) approximately 1000 
cards. 


CIRCUS FORTRAN SOURCE DECK. Approximately 
6000 cards. 


CIRCUS ASSEMBLER LANGUAGE ROUTINES: 
Approximately 500 cards total. 


CIRCUS MODEL LIBRARY UPDATE PROGRAM 
FORTRAN SOURCE DECK, exclusive of those 
routines common to CIRCUS: See CIRCUS PRO- 
GRAM MANUAL (No. 346-2) less then 500 cards. 


CIRCUS SAMPLE DATA. Approximately 500 cards. 


The 14 files fill approximately 1000 feet of tape. 


All files on all tapes are 80 character card images: 80 byte 


records unblocked. 


All tapes are odd parity, density 800 bits/inch. 


Note that the first two files only are in object form, all others 


are in EBCDIC. 
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Ne LES ON GERGUS/ 360: 

In order to CIRCUS/ 360 as supplied, the user must provide 
timing routines which are called by the program. These routines are 
a part of Lincoln Laboratory's system library and therefore are not 
a part of CIRCUS. Dummy routines may be written, but if different. 
names are used, or if similar routines exist in the user's system but 
with different names, the calling statements must be changed in 
CIRCUS. Calls appear in KLOCK and OUTPUT in several places in 
these routines. The routines are: 

TIMDAT (IYR, IDAY) 
TIMHR (HOURS) 
TIMHMS (IHR, IMIN, ISEC) 


where I implies fixed point and 
H implies floating point numbers 


These routines use the system clock and times as implied in their 
arguments. 

If the system on which CIRCUS/ 360 is to be used has a core 
capacity greater than 256 K bytes, the capacity of CIRCUS may be 
increased to handle circuits of greater than 70 equations by adding 

COMMON A1(50000) CIR25715 
and changing CIR25750 to read 


LGGrai. = 50000 CIR25750 


in SUBROUTINE MAPPER. 


70 


In order that the supplied version fit into a 256 K byte core, 
the system's load point must be BOOO or less. 

CIRCUS/ 360 uses the system clock via the assembler language 
routine SETLMT to compare execution time against a specified limit. 
If no limit is specified, the program assumes a ten minute limit, 
after which CHARET is called and execution terminated. To adjust 


this limit, the first card in the data deck should read, beginning in 


column one, 


LIMITS, X 
where X is the desired limit in minutes, floating point (as in LIMITS, 
8.5 for 8 min. 30 sec. maximum execution). 

CIRCUS/ 360 now includes a printer plotter routine which is 
called by the PLOT statement in the input list (see user's manual). 
Plots are 51 by 120 characters; time interval is the smallest in the 
INTERVALS statement; if greater than 120 time points at this interval 
exist, the plot will be continued on succeeding pages. One plot vari- 


able appears per plot. 


KNOWN BUGS IN CIRCUS/ 360 

Be Occasionally a circuit time and associated outputs may be 
printed twice. This occurs at times when something is 
specified in the input list as changing. The first of these 
printed times represents a value slightly less than that 


printed by an amount determined by machine accuracy, and 
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rounded off when printed. The second value is the exact time 
as specified in the input list. 

When specifying a logical unit for a model library, save tape, 
etc., a user is supposed to be able to specify the logical 

unit number (as in SAVE, 15) or use a letter/number 
equivalent (as in SAVE, E) where the letter/number cor- 
respondence is specified in LINK2 (DATA MNSAVE/ 1HA, 
1HB,1HC,]HD,1HE,1HF,1,2,3,4,15,16/ CIR20600). 

As it now stands, the number designation will not be 
recognized, while the letter designation works. If a number 
other than those in the standard letter/number correspondence 
is desired, the aforementioned DATA statement may be 


changed to provide a more useful correspondence. 


(2 


APPEN DIX 


THE CODING OF POTTLE 


LINEAR NETWORK ANALYSIS PROGRAM 


Given an arbitrary, consistent, linear, fixed electrical 
network (with zero initial conditions), the POTTLE program cal- 
culates: 

(a). A prescribed but arbitrary transfer function t(s) or, 

more generally, an arbitrary n-port matrix Cae 

(b). The poles and zeros of t(s) or Py 

(c). The frequency responses associated with ae - 

prints and printer-plots magnitude, phase, and delay. 


(d). A set of state equations for the network. 


In addition, any set of element values may be stepped up to 
six times and the above analysis repeated at each step. 

The network size limitations are: 

(a). Maximum number of elements - 60 

(b). Maximum number of LC elements - 32 

(c). Maximum number of nodes - 27 

(d). Number of controlled sources plus number of inputs 

must be less than or equal to 32. 
(e). Number of controlled sources plus number of outputs 


plus number of resistors must be less than or equal to 64. 
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(f). Three time (number of LC) plus two times (number of 
inputs) plus number of branches minus number of 
nodes must be less than or equal to 128. 


(g). Maximum number of stepped element values - 88 


The input data format is: 


(a). A title card - columns 1-80 will be printed on the 


output sheet. 


(b). Cards describing the branches of the network: 
Columns Des cription 


1-4 Branch name (left justified) beginning 


with: 


V - voltage source 

I - current source 

R, L, or C - resistor, inductor, 
or capacitor, respectfully 


6 - 10 Name of controlling branch (left justified) 
if col 1-4 contains a controlled 
source; otherwise leave blank 

tree lee Positive node (right justified) 
[4 a? (5 Negative node (right justified) 
Column SS 


Column 
ime 2 ee ee Mess 


Number nodes consecutively from zero 
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Columns Description 


18 If branch variable is to be considered 
an output quantity, use symbols V 
or las appropriate. 


20 Type of control (V or JI) if controlled 
source; otherwise blank 


Zi - 30 Value of element or strength of con- 
trolled source; independent source 
leave blank. Must have a decimal 
point (F 10.0 format) 


3a 6 - 6 40 Values of stepped elements or 
strengths of stepped controlled 
42 - %0 sources; all stepped elements 


must be stepped the same number 
of times; stepped element val- 
ues may not be continued on 
72 -. 80 succeeding cards. Must have 
decimal point (F10.0 format) 
(c). A card denoting the types of responses desired. This 


card must be present. 


columns 11 - 12: 01 if frequency response is 
desired 


00 or blank otherwise 


(d). Two cards to describe frequency response desired; these 


are omitted if no frequency response is desired. 


Cava 7]: 
coltmn ¥: 1] if linear frequency scale 
2 if log frequency scale 
column 3: 1 if magnitude calculated in db 


2 if magnitude in ohms 
3 if magnitude in mhos 


©) 


column linear frequency scale 
total number of frequencies 


Oo - 
logrithmic frequency scale 
number of frequencies per 
decade 
number is right justified 
(13 format) 
Caimds2 


columns 1 - 10: Lowest frequency (Hz) 


columns 11 - 20: Highest frequency (Hz) 


(e). If more than one network is to be analyzed, a set of these 


cards is required for each additional network. 


Additional comments: 

(a). The method of analysis used is sensitive to scaling, 
producing overflows if improperly scaled. Values for 
L's and C's should be within an order of magnitude of 
unity for large networks. 

(b). Zero values for inductors and capacitors may be used 
to establish a dummy branch current or node-pair 
voltage; this is useful in controlling a controlled source 


or defining an output variable. 
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APPENDIX D 


THE POTTOL TRANSLATOR 


PURPOSE 

The primary purpose of the Translator is to convert an 
engineer-oriented free-format circuit description into the standard 
input format required by C. Pottle's Linear Network Analysis 
Program. Prior exposure to the POTTLE program (See Appendix 


C) is assumed. 


FEATURES 
The major features of the Translator are: 
(1). Free-format of input data 
(2). Minimal changes to the original program 
(3). Simple and effective overlay (if needed) 
(4). Additional diagnostics 
(5). Exact reprint of input data 
(6). Continuation cards 
(7). Comment cards 


(8). Batch processing 


OVERLAY 
The particular form in which the Translator has been con- 
structed lends itself to a very natural and fast overlay structure. 


Figure D-1 shows the recommended overlay structure. This particular 


(7 





MAIN PROGRAM 
SUBROUTINES 
iN PUT 
DECODE 













TRANSLATOR 


Place all other routines 
in this branch 


SHIFT is in Assembly language 







LANPRO FREPRO 
TNUM OULEA 
PACK TA 






EDIT Sion 





Figure D-1. Overlay Tree Structure. 


structure loads the Translator along with the root segment, processes 
the input data, loads the rest of the program, executes the program 
and stops without having to repetitiously load and reload the same 


overlay branch. 


POTTLE PROGRAM COANGES 

Since the logical unit assignments used in the original program 
were not the same as used locally, it was judged worthwhile to change 
all logical unit references to variables and assign them at one central 
location. This is presently accomplished in the main program of the 
program. These assignments are transferred via labeled COMMON 
to the various subroutines. The present assignments and their 
correlation with the original assignments are shown in figure D-2. 
These assignments may be changed permanently by changing the 


appropriate cards in the main program. Additionally a temporary 
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a 


INTAPE 













User's input data 





IOUTE Copy of input and errors 





STATP State equations 


TRANTP 





Transfer functions, poles/zeros 





PREOTP Frequency response 


Figure D-2. I/O Unit Assignments. 


change of any or all of units may be invoked for a single run through 


the use of a TAPES card as shown in figure D-3. 


column 7 or later 


V 


TAPES, STATP = 2, TAN TP 





| Figure D-3. Temporary Unit Assignments 


ADDITIONS 
The additional features added to the allowable input include: 
(1). Continuation Cards: Any number of continuations are 
allowed as long as the previous card's last character 
is a comma. There is a limit of 154 non-blank charactors 


after column 6 of the first card. 


fa, 





(2). Debugging: A limited debugging facility is provided in» 
that a card containing the word TRACE (starting in 
column 7 or later) will cause a diagnostic print as control 
enters and departs each subroutine. It also provides a - 


print out of the NHEX variables in Al format. 


(3). Batch Processing: Batch processing is implemented 
by requiring that each set of circuit description data 
end with an EXECUTE card and the batch run or job 


end with an END card. 


(4). Comment Cards: Except for reproduction in the print 


out, any card with an asterick in column one is ignored. 


(5). Number Format: All numbers on the input data cards 
may be either F, E, or integer format with or without 
decimals. Recognizable number forms are: 

-XX.X E-YY XX OO: VV 


XXXE+ YY Rex ma ose 


If the decimal is missing it is assumed to follow the 


last X. 


(6). A print out of the translated program is provided. 


INPUT DESCRIPTION 
The standard R,L,C,V, or I branch is coded by placing a 


unique name for the branch in columns 1-5 (note: This is the only 
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type of card allowed to have anything in the first five columns). If 
five characters are used in the name the fifth character will be 
ignored; however, the name need not start in column one. All other 
information is placed after column six (including continued cards). 
The node numbers are delimited by parens and separated by a comma. 
The order of the node numbers implies the positive direction of 
current flow. Spaces may be used freely to enhance readability. 


Typical examples are shown in figure D-4. In the case of stepped 





Figure D-4. Input Data Examples. 


elements the initial value (or only if not stepped) is placed first with 
sequential stepped parameter values separated by commas following. 
All stepped parameters must be stepped the same number of times. 
An output of the branch voltage or current is indicated by adding 
the expression OUTPUT = V or OUTPUT =I, respectfully, to the 
branch description. 

Descriptions of controlled sources have one additional restric- 
tion. The controlling branch name must appear BEFORE the type 


of control. This is illustrated in figure D-5. 
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Figure D-5. Description of Controlled Sources 


FREQUENCY RESPONSE 

The inclusion of a card of the form shown in figure D-6 will 
generate the proper commands to obtain the indicated fA uancy 
calculations. This card will set the variables via labelled common 
required in subroutine TPRINT. It is not necessary to include the 


words LOG and/or DB if these options are desired. 


COMMANDS 
The EXECUTE card causes an analysis to begin. It also 
generates the blank card required by the POTTLE program if no 


frequency card is included. 


FREQ=nl (n2) n3, LINEAR or L@G, @HMS or MH@S or DB 





nl-min freq (Hz Defaults 

n2-no. freq/ decade (log) FREQ none 
no. freq (linear freq scale | LOG 

n3-max freq (Hz) units DB 


Figure D-6. The Frequency Card. 
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An END card terminates the particular run. Lack of an END 
card does not adversely affect the execution of the program (if all 


other data is correct). 


COMMENTS 


The Translator is intimately tied to the POTTLE program. 
Therefore a proper appreciation of the Translator requires at least 
a passing acquaintance with the use of the original POTTLE program. 
To this end Appendix C contains a description of the original program 
and the CADLIB contains the original program. 

In order that POTTLE changes could be kept to an absolute 
minimum the Translator operates on only one card (and its con- 
tinuations) at a time, i.e., the order of input is the same as the 
original POTTLE program except that the label card is no longer used 


as unlimited comment cards are available for expositor comments. 
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APPEMDIRG E 


CALAHAN - A USER'S GUIDE 


CIRCUIT LIMITATIONS 


(1) Complexity 
(a). 30 nodes 
(b). 100 passive elements 
(c). 20 active elements 


(2) Voltage controlled current source is the ONLY 





allowable active element. 


(3) Storage - approximately 90 K bytes 


DEFINITIONS 


(1) Network - the circuit to be analyzed when considered 


asa two-port black box (See figure E-1 below). 






O 
Input | Vi Network V Output 
O 


Figure E-1. The Calahan Network. 


(2) Network function - A functional relationship between the 
input voltage (V,); the output voltage (V,); the input 
current (I); and/or the output current (I,) where these 


terms are positive as indicated in figure E-1. 


84 


(3) KEY 1 - A computer flag whose value indicates a particular 


network function (See Table E-1 below). 


V 
1 Voltage Transfer Function < 
I2=0 

























V 
2 Open Circuit Driving V, 
Point Impedance Sea | 
7 
Open Circuit Transfer V 
—h 
Impedance 
if 
Short Circuit Driving Lr 


Point Impedance 


Short Circuit Transfer 
Admittance 








Current Transfer Function 


Table E-1. Definitions of KEY 1. 


(4) KEY 2 - A computer flag whose value indicates a particular 


set of desired outputs (See table E-2). 


INPUTS 
(1) Network elements or network function (Voltage or current 
transfer function only. 
(2) Input and output node numbers 


(3) Type of network function desired as output 
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KEY 2 | Output(s) is (are) of the Network Function 


Coefficients, poles, and zeros 


Coefficients, poles, Zeros, and freq response 


Coefficients, poles, 2eros, and time response 


Coefficients, poles, zeros, freq response, and 
time response 


Symbolic form of the network function 





te Table E-2. Definitions of KEY 2. 


OUTPUTS (See Table E-2) 


(1) Coefficients of requested network function 

(2) Poles and zeros of the requested network function 

(3) Frequency response (both tables and printer plots of 
magnitude, phase and delay) 

(4) Time response to an arbitrary input (combination step 
in CADLIB version). Tables and printer plots are 
included. 

(5) Repeated outputs of 1-4 with some element value as a 
variable parameter 


(6) Designated network function in symbolic form 


INPUT CEBING 


The initial step in coding for CALAHAN is to number all nodes 


sequentially from one. A deck of cards is then coded as indicated on 
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the flowchart of figure E-2. Two points need clarification. The first 
card is a title card and is ignored except for titling the output. The 
second card is a dual purpose card. If the input is in network 
element form, card #2 contains a positive number equal to number 

of circuit descriptions between title cards. As each circuit or circuit 
variation will normally be different the usual value will be one. If 


the input is to be a transfer function the value used will be a negative 


KEY @. 


DISCUSSION 

General Precautions 

(1). Poles, zeros, and time response may be impossible to 
find for no predictable reason. 

(2). Some bandpass networks may yield incorrect responses 
due to single precision limitations. 

(3). The nodes must be numbered 1! through N, the number 
of nodes; otherwise, there are no constraints on the 
manner of numbering. 

(4). The network must be connected (i.e., a path must exist 
between any two nodes). Thus, for a transformer, the 
common connection between windings must be made, 
unless a path connecting the windings already exists 


through external circuitry. 
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(58 


(6). 


aia) 
(8). 


(9). 


(10). 


(Te 


(12). 


(13) 


The RLC elements must be listed before the controlled 


sources, although in any order. 


Two or more RLC elements across the same nodes must 


have nodes listed in the same order. 

Negative RLC elements are permitted. 

A zero value for a resistance or inductance is never 
permitted. 

If repeated analyses are made with some element value 
as parameter, a zero value is never permitted fora 
resistance or inductance, and is permitted only as a 
final value for a capacitance or controlled source. 
When a logrithmic frequency scale is used, recall that 
zero frequency is not a valid lowest frequency. 

When choosing a highest frequency, beware of possible 
overflow. If a is the coefficient of the highest power 


(n) of s in the denominator, then overflow may occur if 


(2.3 f a yor is greater than (10) ® 


The time response cannot be calculated for a network 
function with a numerator of higher degree than the 
denominator. 

The topological program can be fooled by certain de- 


generate cases which do not warrant safeguarding. 
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